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SECTION  I 


INTRODUCTION 

The  objective  of  this  study  vas  to  Investigate  large  mode  volume 
hole -co up led  optical  resonators  for  high  power  laser  applications.  The 
work  program  was  divided  Into  six  tasks  with  the  following  primary  ob- 
jectives. 

Task  1 — To  Hat-ermine  the  usefulness  of  hole  coupled  resonators  a 
computer  program  was  to  be  developed  to  calculate  and  plot  the  eigen- 
values, the  Internal  resonator  fields,  the  far-fleld  pattern  of  the 
output  beam,  and  the  sensitivity  of  length  misadjustment.  Two  specified 
hole-coupled  resonators  were  to  be  analyzed. 

Task  2 - A computer  program  was  to  be  created  to  optimize  the  figure 
of  laser  resonator  mirrors  such  that  the  dominant  resonator  mode  would 
be  the  best  fit.  In  a least  squares  sense  to  a desired,  user  specified, 
mode.  The  test  case  would  be  a uniform  amplitude  and  uniform  phase 
distribution  over  a circular  aperture  at  the  output  plane  of  the  resonator. 
Once  the  mirror  figures  were  found,  the  resonator  was  to  be  analyzed 
using  the  tools  developed  In  Task  1. 

Task  3 - A study  was  to  be  made  to  determine  the  ratio  of  annular 
to  hole  output  coupling  required  to  achieve  good  azimuthal  mode  discrim- 
ination In  a designated  resonator  configuration. 

Task  4 - A study  similar  to  that  of  Task  3 was  to  be  conducted  for  a 
resonator  with  holes  in  both  mirrors. 

Task  5 - An  investigation  was  to  be  made  of  the  possibility  of  using 
conformal  mapping  techniques  to  transform  circular  geometries  to  rectan- 
gular geometries  where  the  resonator  eigenmode  calculations  could  be 
expedited  using  the  fast  Fourier  transform. 


Task  6 - The  field  fitting  program  of  Task  2 was  to  be  improved 
and  applied  to  various  "desired''  mode  shapes  while  noting  the  effects 


on  the  azimuthal  mode  discrimination  properties  of  the  associated  resonators. 

The  results  of  Tasks  1 through  5 are  presented  in  Sections  II 
through  VI,  respectively.  Since  the  results  of  Task  6 primarily  yielded 
improvements  in  the  results  of  Task  2,  they  are  included  in  Section  III. 

The  capabilities  of  two  computer  programs  developed  during  the 
contract  period  are  described  in  Sections  II  and  III.  Details  related 
to  the  calculation  methods  used  in  the  programs  are  contained  in 
appendixes  A and  B. 
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SECTION  II 


TASK  I - RESULTS 

A.  General  Information 

A computer  program  has  been  developed  and  used  to  analyse  the  eigen- 
modes  of  several  hole-coupled  resonators.  The  program  can  compute  and 
plot:  1)  the  eigenvalues  versus  the  outer  resonator  Fresnel  number,  2) 
the  internal  resonator  fields  (magnitude  and  phase),  and  3)  the  far-field 
patterns  of  circular  cylinderical  resonators  having  mirrors  of  arbitrary 
radial  figure  in  the  Fresnel  approximation. 

The  method  of  calculation  used  in  the  program  is  similar  to  that 
described  by  Siegman  and  Miller  in  their  article  on  the  Prony  method  [1]*, 
but  with  appropriate  modifications  to  accommodate  arbitrary  radial  mirror 
figures.  To  dispense  with  the  problem  of  the  "spurious"  eigenvalues 
which  they  described,  the  Prony  calculation  itself  is  iterated  to  insure 
converged  (or  non-converged)  solutions.  The  initial  vector  for  each 
Prony  calculation  after  the  first  one  is  taken  as  the  sum  of  the  nor- 
malized eigenvectors  computed  during  the  previous  calculation.  In  cases 
where  convergence  is  not  achieved  after  a specified  number  of  iterations, 
those  eigenvalues  and  eigenfunctions  are  discarded.  More  specific  details 
concerning  the  method  of  calculation  and  computer  programming  are  contained 
in  Appendix  A. 

Two  hole-coupled  resonators  were  identified  by  the  Air  Force  project 
officer  [2]  as  test  cases  for  analysis  using  the  computer  program.  The 
mirror  profiles  of  these  two  resonators,  designated  HUR  and  HUR  CC  are 
depicted  in  Figure  1.  The  other  two  resonators  shown  in  the  same  figure 
are  considered  "negative  branch"  versions  of  the  first  two.  These  four 
resonators  are  confocal  versions  of  a more  general  class  of  resonators 


^Numbers  in  brackets  indicate  references. 


In  which  the  mirrors  are  conic  sections  with  arbitrary  off-axis  spherical 
curvature.  These  are  "confocal"  but  not  in  the  usual  sense  because  the 
mirrors  have  focal  rings  rather  than  focal  points.  The  radii  of  the 
focal  rings  are  designated  and  D2  for  mirrors  1 and  2,  respectively. 

Both  mirrors  have  an  outer  radius  of  A.  Mirror  2 has  a central  output 
aperture,  or  hole,  of  radius  H.  A magnification  M ” A/(A-H)  is  asso- 
ciated with  the  first  two  resonators,  and  a negative  magnification 
M ■ A/ (H-A)  is  associated  with  the  "negative  branch"  configurations. 

The  computer  program  was  used  to  plot  the  radial  mode  eigenvalues 
versus  the  outer  resonator  Fresnel  number  N • A^/(BX)  (B-mirror  separa- 
tion, X-wavelength)  for  various  magnifications  and  azimuthal  mode  indices 

1 (where  azimuthal  variation  of  the  form  exp  (+  j£<J>)  has  been  assumed). 

The  program  was  also  used  to  plot  the  internal  resonator  fields  and  the 
far-field  patterns  for  various  Fresnel  numbers.  Representative  plots 
are  shown  and  discussed  in  the  following  subsections. 

B.  Analysis  of  Resonator  1 - HUR 

Figures  2 through  6 show  the  eigenvalue  plots  for  resonator  1 of 
Figure  1 for  magnifications  M - /l,  5/3,  2,  3,  and  8,  respectively.  The 
eigenvalues  intertwine  in  a slightly  more  complicated  fashion  than  those 
of  the  conventional  unstable  resonator  (e.g.,  see  reference  1).  An  even 
more  significant  difference  between  the  mode  of  this  resonator  and  a conven- 
tional unstable  resonator  is  that  only  for  very  small  Fresnel  numbers 
(N-l)  is  the  1*0  mode  dominant.  (The  dominant  mode  at  a given  Fresnel 
number  is  the  mode  with  the  largest  eigenvalue  magnitude  |y(  or  lowest 
losses  1 - |y J 2 • ) For  example,  in  comparing  parts  a,  b,  c,  and  d of  Figure 

2 for  M ■ /l,  the  Zm0,  1,  2,  and  3 azimuthal  modes  are  dominant  within 
the  approximate  Fresnel  number  ranges  0-1.5,  1.5-5. 5,  5. 5-9. 5,  and  9.5-?, 
respectively. 
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Figure  6d  A ■ 3 


Figure  7 shows  the  radial  eigenmode  of  field  plots  for  the  lowest 
order  azimuthal  mode  lm0 , for  M**2  and  N«6.  The  eigenvalues  for  the  three 
radial  modes  can  be  found  In  Figure  4a  at  N-6.  Figure  7a  shows  the 
amplitude  and  phase  of  the  three  lowest  loss  radial  modes  In  a plane 
just  after  being  reflected  from  mirror  1.  Probably  the  most  striking 
feature  Is  the  rather  sharp  magnitude  peak  at  the  center  of  the  mirror 
(relative  radial  coordinate  equal  to  zero).  The  phase  of  the  lowest  loss 
radial  mode  decreases  approximately  90°  from  the  center  of  the  mirror  to 
the  outer  edge  (relative  radial  coordinate  equal  unity).  The  phase  of  the 
next  lowest  loss  mode  on  the  other  hand  Increases  to  approximately  90° 
near  the  outer  edge  but  then  decreases  to  near  zero  at  the  edge.  The 
phase  profile  of  the  third  lowest  loss  mode  has  even  more  variation. 

Figure  7b  shows  the  same  fields  except  plotted  In  a plane  just 
Incident  on  mirror  2.  At  this  magnification  the  radius  of  the  hole  In 
mirror  2 Is  one  half  the  outer  mirror  radius.  Since  the  fraction  of 
the  field  In  the  hole  relative  to  that  on  the  mirror  Increases  for  the 
higher  order  radial  modes  It  Is  easy  to  understand  that  the  higher  order 
modes  have  more  losses.  Both  the  magnitude  and  the  phase  of  the  lowest 
loss  radial  mode  are  fairly  uniform  in  the  output  hole  which  results 
In  the  near  "ideal"  far-field  beam  profile  shown  In  Figure  7c.  For 
comparison,  the  far-field  pattern  corresponding  to  uniform  illumination 
of  the  hole  with  the  same  power  as  in  the  lowest  order  radial  mode  for 
1*0  is  also  shown , Finally  the  integrated  far-field  Intensity  associated 
with  the  two  far-field  plots  of  Figure  7c  are  shown  In  Figure  7d.  Out  to 
the  half-power  beam  width,  the  two  plots  are  virtually  identical.  The 
fact  that  there  Is  slightly  more  power  within  a given  beam  angle  larger 


than  the  beam  width  for  the  non-uniformally  illuminated  aperture  relative 
to  the  unlformally  illuminated  aperture  is  probably  not  too  significant, 
since  the  power  is  reduced  there. 

Figures  8,9,  and  10  show  the  fields  for  the  same  conditions  as 
Figure  7 for  the  azimuthal  mode  Indices  1 ■ 1,  2,  and  3.  For  these 
higher  order  azimuthal  modes,  the  field  is  zero  on  the  axis.  Also,  it 
can  be  seen  that  the  fields  tend  to  "move"  radially  outward  for  increasing 
azimuthal  mode  index.  In  comparing  the  field  magnitude  profiles  at 
mirror  2 for  the  1* 0 and  1*3  dominant  radial  modes  (Figure  7b  and  10b) , 
it  is  understandable  that  the  1*3  mode  has  less  loss  than  the  1*0 
mode  since  there  is  relatively  little  field  in  the  hole  for  the  1*3  mode. 

For  all  the  magnifications  considered,  the  1*0  modes  were  dominant 
only  at  very  small  Fresnel  numbers.  Increasing  the  size  of  the  hole 
from  zero  apparently  increases  the  losses  more  rapidly  for  the  1*0  mode 
than  the  higher  order  azimuthal  modes.  Thus,  to  obtain  any  significant 
power  output  from  this  type  resonator  in  the  1*0  mode  will  require  that 
some  method  be  devised  which  would  discriminate  against  the  higher  order 
azimuthal  modes. 

C.  Analysis  of  Resonator  2 - HUR  CC 

Figure  11  shows  the  eigenvalues  for  the  four  lowest  order  azimuthal 
modes  for  resonator  2 of  Figure  1 with  M*2.  Only  for  Fresnel  numbers  less 
than  0.5  is  the  1*0  mode  dominant.  The  computer  program  has  considerable 
difficulties  in  achieving  converged  solutions  for  Fresnel  numbers  greater 
than  N*6.  A clue  to  the  difficulty  can  be  obtained  by  examining  the 
complex  field  structure  for  the  1*0  and  1*1  modes  in  Figures  12  and  13 
at  N*4.  The  large  number  of  grid  points  required  to  delineate  these 
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Figure  8b  Radial  mode  profiles  in  a plane  incident  at  mirror  2 
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Figure  9b  Radial  mode  profiles  in  a plane  lncldenc  at  mirror  2 
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Far-field  pattern  of  lowest  loss  radial  moae 
and  reference  field  (uniformly  illuminated 
aperture) 
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Figure  12a  Radial  mode  profiles  In  a plane  reflected  from 
mirror  1 of  HUR  CC  for  M ■ 2,  A ■ 0,  N ■ 4 
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"interference"  patterns  makes  the  calculations  at  Fresnel  numbers  larger 
than  about  N«6  prohibitively  expensive.  Even  if  the  )l*0  were  dominant 
at  reasonably  large  Fresnel  numbers,  this  resonator  would  be  objectionable 
because  of  the  poor  far-field  performance  as  seen  in  Figure  12d.  From 
these  and  other  calculations  it  has  been  concluded  that  this  resonator 
configuration  would  not  provide  a high  quality  output  beam. 


D.  Analysis  of  Resonators  3 and  4 

Even  though  resonators  3 and  4 were  not  specifically  designated 
for  analysis,  some  preliminary  calculations  were  made  to  determine  their 
potential  for  use  as  hole-coupled  resonators.  Figures  14,  15,  and  16 
show  the  eigenvalues  calculated  for  resonator  3 of  Figure  1 for  magni- 
fications  M-  -1,  -2,  and  -3.  These  plots  are  interesting  in  that  the 
radial  mode  eigenvalues  do  not  cross  as  has  been  the  case  for  all  pre- 

% 

viously  examined  stable  and  unstable  resonators.  Another  interesting 
and  potentially  more  important  aspect  of  these  graphs  is  that  the  azimuthal 
mode  eigenvalues  tend  to  become  degenerate  as  the  Fresnel  number  increases. 

Thus  only  minimal  azimuthal  mode  discrimination  should  be  required  to 
achieve  £*0  mode  operation.  Hence  if  the  theory  of  reference  [3]  about 
the  role  of  mode  volume  with  regard  to  mode  discrimination  by  an  active 
medium  is  correct,  then  the  2*0  mode  could  be  the  dominant  mode  when  the 
resonator  contains  a saturable  gain  medium  since  the  2*0  mode  is  the  only 
mode  that  does  not  have  zero  amplitude  on  axis. 

The  excellent  radial  mode  discrimination  and  the  possibility  that 
2*0  mode  may  be  favored  because  of  its  mode  volume  make  resonator  3 the 

I 

most  promising  hole-coupled  resonator  investigated  to  date.  In  addition 

! 

to  continued  theoretical  analysis,  an  experimental  study  of  this  type  of 

1 

resonator  should  be  conducted. 
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A few  calculations  of  the  eigenvalues  of  resonator  4 (not  shown) 
revealed  that  the  higher  order  azimuthal  modes  have  lower  losses  than 
the  £-0  mode.  Therefore  no  further  analysis  of  this  configuration  seems 
worthwhile,  since  we  are  primarily  interested  in  only  those  resonators 
where  the  4*0  mode  would  be  dominant. 


E.  Sensitivity  to  Length  Misadjustment 

One  method  to  determine  the  sensitivity  of  the  eigenmode  character- 
istics to  length  misadjustment  would  be  to  calculate  the  eigenmodes  for 
different  distances  between  the  mirrors  with  all  other  parameters  fixed. 

A second  method  is  to  note  that  changing  the  resonator  length  a small 
amount,  to  first  order,  just  changes  the  resonator  Fresnel  number.  Hence, 
to  determine  for  example  the  sensitivity  of  the  eigenvalue  magnitude 
| Y | to  mirror  separation  B,  one  would  compute 
S I Y I _ B | IyI  a B AixL 

B TyT  3 b TyT  AB 

where  the  A signifies  a smell  change  in  the  respective  parameter.  Note 

Ail!  , All!  AS.  . S Ail! 

AB  AN  AB  B AN 

since  the  Fresnel  number  N is  Inversely  proportional  to  B.  Thus 

s IyI  « n Aii!  . . ~ lYi‘ 

B TyT  an  Ty^T  n2  - Nx 

where  (from  the  eigenvalue  graphs)  |y2I  and  |y2I  are  the  eigenvalue 
magnitudes  at  Fresnel  numbers  N^  and  N2  corresponding  to  the  mirror 
separations  B and  B + AB,  respectively.  In  a similar  manner  other  sen- 
sitivities related  to  Fresnel  number  changes  can  be  computed. 
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SECTION  III 


RESULTS  OF  TASKS  2 AND  6 

A.  Introduction 

A computer  program  has  been  developed  which  calculates  the  resonator 
mirror  figures  to  obtain  an  optimum  fit  in  a least  squares  sense  to  a 
specified,  or  desired,  radial  eigenmode  amplitude  and  phase  profile. 

A description  of  the  calculation  technique  used  in  the  computer  program 
is  given  in  Appendix  B.  The  program  was  to  be  tested  by  noting  how  well 
it  performed  in  trying  to  achieve:  1)  a maximum  mode  volume,  2)  uniform 
amplitude  and  phase  in  the  output  of  a hole-coupled  resonator,  and  3) 
maximum  far-field  on-axis  intensity.  Tests  conducted  during  work  on 
Task  2 indicated  that  the  program  performed  satisfactorily  except  when 
the  resonator  had  little  mode  discrimination.  Later  while  working  on 
Task  6,  the  Prony  method  of  calculating  the  fields,  which  was  used  in 
Task  1,  was  incorporated  into  the  program,  and  the  program's  performance 
improved  significantly.  The  test  cases  were  repeated,  and  better  fits 
to  the  desired  fields  were  obtained.  These  latter  tests  results  are 
described  in  the  following  subsections,  along  with  some  additional  test 
results  obtained  during  Task  6. 

Before  results  are  described,  the  main  limitation  of  this  approach 
to  designing  resonators  should  be  mentioned.  The  program  calculates  the 
mirror  curvatures,  or  figure,  for  an  optimum  fit  of  the  lowest  loss  radial 
mode  of  a particular  azimuthal  mode  index  to  the  desired  mode  profile. 

In  all  the  test  cases,  the  i* 0 mode  index  was  used.  After  the  "optimum" 
resonator  had  been  found,  it  was  analyzed  using  the  tools  of  Task  1.  In 
almost  all  cases,  a higher  order  azimuthal  mode  was  found  to  have  fewer 
losses  than  the  im0  mode  for  which  the  resonator  was  designed.  To 
assure  that  the  resonator  is  designed  for  the  lowest  loss  radial  and 
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azimuthal  mode,  the  one  dimensional  radial  field  propagator  would 
have  to  be  replaced  with  a full  mirror  propagator.  Neither  time 
nor  sufficient  computer  funds  were  available  to  make  this  important 
modification. 

B.  Maximum  Mode  Volume 

In  an  attempt  to  achieve  a maximum  mode  volume,  the  program  was 
given  the  task  of  trying  to  achieve  a uniform  field  magnitude  over  the 
entire  output  plane  of  the  hole-coupled  resonator  shown  in  Figure  17. 
Maximizing  the  field  in  one  "plane"  does  not  guarantee  a maximum  mode 
"volume",  but  the  program  is  currently  structured  to  fit  fields  in  only 
one  plane.  Since  the  least  squares  minimization  procedure  finds  a 
local  minimum  for  the  mean  square  error,  it  is  best  that  the  initial 
resonator  configuration  be  one  in  which  the  fields  are  already  close  to 

j 

the  desired  fields. 

The  beginning  and  ending  resonator  configurations  along  with  the 
eigenvalues  for  the  two  lowest  order  radial  and  four  lowest  order  azimuthal 
modes  indices  (m-0,1  and  1-0, 1,2, 3)  are  shown  in  Figure  17a.  Figure  17b 
shows  the  desired,  beginning,  and  ending  field  on  the  output  mirror 

(2)  for  1-0.  Note  that  the  ending  field  magnitude  is  a very  good  fit 

1 , , 

to  the  desired  field  but  that  the  eigenvalue  for  this  mode  ( |Yqq!"0.656) 

j is  not  the  largest  eigenvalue.  For  reference  purposes  the  associated 

fields  on  the  other  mirror  (1)  are  shown  in  Figure  17c. 

It  appears  that  the  ending  field  on  mirror  1 has  less  mode  "volume" 
than  the  beginning  field,  but  the  situation  is  somewhat  distorted  by  the 
central  peak  since  the  fields  are  normalized  by  peak  amplitudes  rather 
than  average  intensity  as  would  be  required  for  a fair  mode  volume  com- 
parison. Apparently  the  uniform  fill  on  mirror  2 was  obtained  somewhat 


Figure  17a  Schematic  and  eigenvalues  of  resonator  to  achieve 
a uniform  field  magnitude  on  mirror  2,  N ■ 4,  Z » 
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Figure  17c  Fields  In  s plane  reflected  from  mirror  1 
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at  the  expense  of  the  fill  on  mirror  1.  The  far-f leld  beam  patterns  of 
Figure  17d  were  calculated  for  the  fields  In  the  output  aperture  normalized 
to  the  same  power  level  In  the  aperture.  The  relative  peak  values  given 
are  measured  relative  to  the  uniformly  Illuminated  (uniform  amplitude  and 
phase)  aperture  field.  The  plot  symbol  + for  the  uniformly  Illuminated  case 
Is  somewhat  obscured  by  the  plot  symbol  0 used  for  the  Initial  field. 

In  any  event,  the  far-f leld  pattern  of  the  beginning  resonator  is  better 
than  that  of  the  ending  resonator. 

In  an  attempt  to  improve  the  far-fleld  pattern  while  maintaining  a 
maximum  mode  "volume",  an  additional  constraint  was  Imposed  on  the  desired 
field:  that  its  phase  be  uniform  in  the  output  aperture.  The  beginning 
(same  as  previous  case)  and  ending  resonator  mirror  configurations  are 
shown  in  Figure  18a.  The  desired,  beginning,  and  ending  fields  on  mirror 
2 are  shown  in  Figure  18b.  The  final  field  magnitude,  as  might  be 
expected,  is  not  as  good  a fit  to  the  desired  magnitude  due  to  the 
additional  constraint  on  the  phase  of  the  field  in  the  aperture.  The 
phase  fit  in  the  aperture  is  good  and  results  in  an  almost  ideal  far-field 
pattern  as  shown  in  Figure  18d.  Again,  for  reference  purposes  the  be- 
ginning and  ending  fields  on  mirror  1 are  shown  (Figure  18c).  The 
ending  field  "fills"  the  mirror  very  well  except  for  the  small  region 
near  the  axis. 

C.  Uniform  Magnitude  and  Phase  in  the  Output  Hole 

A test  of  the  computer  program  was  made  in  an  attempt  to  achieve 
uniform  magnitude  and  phase  in  the  output  aperture  of  a hole-coupled 
resonator.  The  initial  resonator  configuration  was  chosen  as  that  of 
the  HUR  described  in  Task  1 (resonator  1 of  Figure  1).  The  beginning 
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and  ending  mirror  profiles  are  depicted  In  Figure  19a.  Note  that  even 
though  mirror  2 was  Initially  convex  with  respect  to  the  Interior  of 
the  resonator,  it  ended  as  concave.  The  fields  for  the  beginning  and 
ending  resonators  at  mirror  1,  at  mirror  2,  and  in  the  far-field  are 
shown  in  plots  b,c  and  d.  As  can  be  seen,  the  fit  to  the  desired  field 
on  mirror  2 is  excellent,  the  mode  fill  is  increased  over  that  of  the 
HUR  on  both  mirrors,  and  the  far-field  is  virtually  identical  with 
that  of  the  uniformly  illuminated  aperture  case.  The  associated  in- 
tegrated far-field  intensities,  "power  in  the  bucket"  curves  are 
shown  in  plot  19 e.  The  total  power  within  the  half-power  beam  width 
for  the  two  cases  are  about  equal;  however,  for  all  beam  angles  greater 
than  the  half— power  beam  angle,  the  HUR  integrated  intensity  exceeds 
that  of  the  uniformly  illuminated  aperture.  The  "tapered"  magnitude 
profile  reduces  the  power  in  the  side-lobes  relative  to  the  main  lobe. 

Next,  a test  was  made  to  sec  if  good  results  could  be  achieved  for 
a resonator  with  a Fresnel  number  of  IM  which  was  twice  that  used  for 
the  previous  tests.  The  beginning  and  ending  resonator  conflgurationa 
and  the  associated  field  plots  are  shown  in  Figures  20a  through  20e.  As 
can  be  seen  the  results  are  very  good.  For  even  larger  Fresnel  number 
calculations  it  would  probably  be  necessary  to  let  the  program  calculate 
the  mirror  curvature  to  higher  (than  second)  order. 

D.  Maximum  On  Axis  Far-Fleld  Intensity 

From  an  analysis  of  the  far-field  Integral  equation  for  a circular 
aperture  it  can  be  shown  that  for  given  power  flow  through  a fixed  aperture 
the  mny-tmum  on-axis  far-field  intensity  is  achieved  when  the  field  in  the 
aperture  is  uniform.  Therefore,  it  was  decided  to  have  the  program  attempt 
to  maximize  the  power  within  a given  beam  angle,  rather  than  to  maximize  the 
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Figure  19a  Schematic  and  eigenvalues  of  a resonator  to 
achieve  uniform  magnitude  and  phase  in  hole 
of  mirror  2,  N - 4,  Si  ■ 0 
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Figure  19b  Fields  In  a plane  Incident  on  mirror  2 
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Figure  19e  Integrated  far-field  intensities 
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Figure  20b  Fields  in  a plane  incident  on  mirror  2 
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Figure  20c  Fields  in  a plane  reflected  from  mirror  1 
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far-field  on-axis  Intensity.  The  approach  taken  was  to  request  that  the 
program  find  a resonator  whose  integrated  far-field  intensity  was  equal 
to  its  maximum  possible  value  over  a specified  beam  angle.  Of  course 
the  ideal  far-field  magnitude  profile  which  would  satisfy  this  condition 
would  be  an  on-axis  delta  function.  The  beginning  and  ending  resonator 
configurations  are  shown  in  Figure  21a.  Figure  21b  shows  the  resulting 
integrated  far-field  intensity,  and  it  can  be  seen  that  the  final  intensity 
is  a better  fit  to  the  desired  intensity  than  that  of  the  uniformly  ill- 
uminated aperture.  This  results  in  a significant  decrease  in  the  first 
side  lobe  peak  relative  to  the  reference  as  shown  in  Figure  21c.  In 
Figure  21d  it  can  be  seen  that  the  field  magnitude  in  the  aperture  of 
mirror  2 is  tapered  and  the  phase  is  uniform.  The  final  field  on  mirror 
1 in  Figure  21e,  although  Itself  quite  peaked  on  axis,  is  more  slowly 
varying  than  the  fields  of  the  other  hole— coupled  configurations  examined. 


E.  A'  "Hole-Coupled"  Resonator  with  a Dominant  2-0  Mode 

Many  additional  tests  have  been  made  of  the  program's  ability 
to  fit  desired  fields,  and  in  most  tests  the  program  was  able  to  achieve 
a good  fit.  Once  the  resonator  parameters  were  found,  the  eigenvalues  for 
the  higher  order  azimuthal  modes  were  calculated  using  the  program  developed 
during  Task  1.  It  was  found  that  none  of  the  hole-coupled  resonators  had 
good  azimuthal  mode  discrimination.  The  only  hole— coupled  resonators 
found  for  which  the  2“0  mode  was  dominant  had  very  high  losses  over  the 
outer  edges  of  the  mirror (s). 

One  "hole-coupled"  resonator  in  which  the  2« 0 mode  was  dominant 
is  shown  is  Figure  22.  The  mirror  profiles  are  extremely  exaggerated 
to  show  their  general  shape.  No  R and  D values  are  given  because  the 
m2rror  figures  were  calculated  to  third  order.  The  beginning  configuration 
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Figure  21a  Schematic  and  eigenvalues  of  resonator  to 
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Is  actually  tha  ending  configuration  of  a previous  test  where  the 
desired  field  was  specified  as  uniform  phase  in  a plane  incident 
on  mirror  2 out  to  r«  ^ where  is  the  radius  of  mirror  1.  No 
restriction  was  placed  on  the  field  magnitude.  In  this  case  the  re— 
quir«Mnt  that  the  phase  be  uniform  on  the  mirror  was  removed.  Except 
near  the  edge  of  the  hole  where  the  field  amplitude  is  relatively  small, 
the  phase  of  the  final  field  is  very  close  to  the  desired  field  phase. 


The  mode  discrimination  ratio 


|I00|  0.515 

'YjO1  0.478 


1.08 


was  the  largest  obtained  for  any  hole— coupled  resonator  examined,  except 
where  the  Fresnel  number  of  the  hole  was  trivially  small.  Even  in  t-M» 
case  the  Fresnel  number  of  the  hole  is  only 

NH1  - (1/3) 2N  - 1 

and  less  than  half  of  the  total  output  is  through  the  hole. 


F.  Rimmed  Resonator  Test 

One  additional  test  will  be  described  to  further  illustrate  the 
utility  of  the  program.  The  program  was  given  the  task  of  adjusting 
the  slope  and  curvature  of  the  rim  on  one  of  the  resonators  described 
by  Lax  et  al.  [3].  They  state  that  "...transverse  mode  can  be  suppressed 
at  high  power  levels,  provided  the  ratio  of  the  fundamental-mode-lrradlance 


spatial  minimum  to  maximum  exceeds  a 


value".  Using  the  fitting 


program  the  fundamental  mode  lrradlance  can  be  tailored  very  easily. 

The  beginning  configuration  (Figure  23a)  was  chosen  without  regard  to 
their  suggestion  that  the  mirror  and  rim  Fresnel  numbers  be  non-integral 
or  that  Nj  ~ N.  We  note  that  there  is  little  azimuthal  discrimination 


since 


iY00|  _ 0.985 
Yin  0.974 
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for  Che  starting  resonator.  The  field  profile  (Figure  23b)  of  the 
4-0  node  for  the  beginning  resonator  Is  interesting  In  that  the  phase 
In  the  plane  at  the  mirror (s)  is  almost  perfectly  flat.  The  desired 


!*< 


magnitude  was  chosen  to  try  to  make  the  field  amplitude  more  uniform 

over  the  mirror  without  increasing  the  losses  too  much,  and  the  desired 

phase  is  that  of  a spherical  wave  coming  from  the  virtual  focus  (see 

reference  4)  of  the  geometric  resonator  mode  with  no  rim.  The  final  field 

magnitude  profile  is  "flatter,"  and  the  aslmuthal  mode  discrimination  is 

somewhat  improved  in  the  ending  resonator,  but  at  the  expense  of  higher 

losses  over  the  mirror  edges.  It  appears  that  the  magnitude  fit  is  not 

as  good  as  might  be  expected,  probably  due  to  the  non-uniform  magnitude 

and  phase  weighting  functions  of  W and  W . The .magnitude  fit  could  be 

m p 

improved  at  the  expense  of  the  phase  fit  simply  by  increasing  the  W weights 

B 

relative  to  the  Wp  weights.  It  is  not  understood,  however,  why  the  phase 
profile  is  not  an  even  better  fit  than  that  shown.  It  is  interesting 
that  the  rim  was  virtually  straightened  out  to  conform  to  the  curvature 
of  the  central  portion  of  the  resonator.  In  any  event,  the  ratio  of 
the  minimum  to  the  maximum  field  values  is  0.3  and  hence  satisfies  the 
criteria  stated  in  [3]  for  achieving  good  azimuthal  mode  discrimination. 

G.  Conclusions 

In  conclusion,  the  program  is  able  to  calculate  mirror  curvatures 
to  fit  a desired  radial  mode  profile  for  no  azimuthal  variations.  It 
has  been  found  that  requiring  the  field  magnitude  of  the  4-0  mode  to  be 
uniform  in  the  output  aperture  virtually  assures  that  the  4- 0 mode  will 
not  be  the  lowest  loss  mode  (l.e.  that  higher  order  azimuthal  modes  will 
have  lower  losses  in  the  same  resonator).  This  suggests  that  possibly 
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Figure  23b  Fields  In  a plans  at  mirror  1 (or  2).  Wm  is  the 
magnitude  error  weighting  function  and  Wp  la  the 
phaae  error  weighting  function  (weighting  function 
described  In  Appendix  B) 
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the  prograa  should  be  required  to  achieve  uniform  amplitude  In 
the  hole  for  the  lowest  loss  azimuthal  mode  to  Increase  Its  losses 
relative  to  the  AH)  mode.  The  "brute  force"  way  to  obtain  A-0 
mode  operation  would  be  to  have  the  program  adjust  the  azimuthal  mirror 
profile  In  addition  to  the  radial  profile.  This  procedure,  while  straight- 
forward with  regard  to  generalizing  the  prograa  to  calculate  the  full 
resonator  field,  would  be  more  costly  with  regard  to  the  computer  time 
required. 


SECTION  IV 


TASK  3 RESULTS 


The  objective  of  Task  3 was  to  detenalne  the  required  ratio  of 
annular  to  hole  coupling  such  that  the  resonator  of  Figure  24  would 
have  good  azimuthal  mode  discrimination.  To  facilitate  the  discussion, 
we  define  the  geometric  ratio  of  the  area  of  the  annular  output  region 
irA^Z  _ ita^  to  the  area  of  the  hole  irH^ 


R 


A22-Al2 


In  the  limiting  case  where  focal  ring  radius  D-0  the  resonator  degenerates 
Into  a conventional  unstable  positive  branch  confocal  resonator  which  is 
known  to  have  good  azimuthal  mode  discrimination.  In  this  case  H^-O  and 
R - In  another  limiting  case  where  D-A^  the  resonator  degenerates  into 
resonator  1 of  Figure  1,  which  was  Investigated  In  Task  1 and  found  to 
have  poor  azimuthal  mode  discrimination.  In  this  second  case  the  area 
of  the  annulus  Is  zero  so  that  R-0.  As  the  area  of  the  hole  increases 
and  R decreases  the  azimuthal  mode  discrimination  ratio  (magnitude  of 
the  ratio  of  the  largest  it— 0 eigenvalue  to  the  next  largest  eigenvalue 
for  any  l)  will  decrease.  Somewhere  within  the  range  0 < R < 00  the  loss 
for  the  lm0  mode  will  equal  the  loss  of  one  or  more  of  the  higher  order 
azimuthal  modes,  and  the  mode  discrimination  ratio  will  be  unity.  It 
was  hoped  that  a large  fraction  of  the  output  could  be  hole  coupled  (R  « 1) 
before  this  breakeven  point  was  reached.  To  determine  how  much  hole  coupling 
could  be  allowed,  while  maintaining  acceptable  azimuthal  mode  discrimination 
the  computer  program  developed  under  Task  1 was  used  to  calculate  the 
eigenvalues  of  the  resonator  for  various  values  of  D. 
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Figures  25  and  26  show  the  eigenvalues  for  the  case  where  D-0 
for  magnifications  M « /?  and  3.  As  Is  seen,  at  a fixed  magnification, 
the  1-0  mode  has  the  lowest  losses  (the  largest  eigenvalue),  and  the 
losses  Increase  as  1 Increases.  It  Is  evident  In  comparing  Figure  25a 
and  25b  that  the  mode  discrimination  ratio  decreases  as  the  Fresnel 
number  Increases. 

Figure  27  shows  the  lowest  loss  eigenvalues  for  D approximately 
one  third  of  the  radius  of  mirror  1 and  for  a magnification  M-3.0. 

In  this  case  the  radius  of  the  hole  ” 2/9A^  and  the  Inner  radius  of 
the  output  annulus  — (5/9)  Al.  Hence 

1 - X5/9?2 

R " ((1/3)2)  " 14 

which  means  that  only  about  l/14th  the  total  output  could  be  hole-coupled. 
A comparison  of  the  different  azimuthal  mode  plots  shows  that  the 
1-0  mode  Is  dominant  In  the  Fresnel  number  Interval  0 < N <4,  the 
1-1  mode  Is  dominant  for  4 < N <15,  and  the  1-2  mode  is  dominant  15  < N < 
a higher  fresnel  number  than  Is  on  the  plot. 

The  eigenvalues  were  also  calculated  for  D - Al/2  for  magnifications 
M - /T  and  3,  as  shown  In  Figures  28  and  29.  As  in  the  previous  case 
the  azimuthal  Index  of  the  dominant  mode  increases  as  the  Fresnel  number 
increases.  Thus  it  is  concluded  that  only  for  relatively  small  Fresnel 
numbers  could  any  of  the  resonator  output  be  hole-coupled  in  the  1-0 
mode.  These  Fresnel  numbers  are  too  small  to  be  of  interest  for  high 
power  laser  applications. 
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SECTION  V 


i 


TASK  4 RESULTS 

The  purpose  of  Task  4 was  to  study  the  behavior  of  the  eigeumodes 

of  the  resonator  shown  in  Figure  30.  This  resonator  is  similar  to  the 

one  in  Figure  24  which  was  analyzed  in  Task  3 except  the  centerline 

has  in  effect  been  "blown  up"  into  a circular  cylinder  of  radius  HI. 

This  resonator  could  be  used  when  there  is  an  on-axis  obscuration 

within  the  region  r A1  caused,  for  example,  by  the  nozzle (s)  of  a 

chemical  laser.  For  parameterization  purposes  the  Fresnel  numbers 

NH1  - H1^/(BA)  and  ND  ■ D^/(BX)  were  defined.  The  behavior  of  the 

eigenmode8  was  analyzed  for  various  values  of  NH1  and  ND  as  a function  of 

2 

the  outer  resonator  Fresnel  number  N • A1  /(BA). 

Figure  31  shows  the  resonator  eigenvalues  for  ND  ■ 1,  NH  ■ 1,  and 
M ■ 2.  This  choice  of  parameters  results  in  a conventional  unstable 
confocal  resonator  whose  centerline  has  been  blown  up  from  a radius 
of  zero  to  a radius  of  HI  which  corresponds  to  a Fresnel  number  of  1. 

A a seen  the  dominant  radial  modes  for  2*0  and  2*1  have  essentially 
identical  eigenvalues  over  the  range  of  Fresnel  numbers  considered 
2 < N < 10.  The  losses  of  the  1*2  and  1*3  modes  Increase  relative 
to  those  of  the  2*0  mode  for  increasing  N. 

Figure  32  shows  the  resonator  eigenvalues  for  ND  and  NH1  increased 
to  a value  of  2 for  M ■ 2.  This  increase  in  HI  results  in  less  azimuthal 
mode  discrimination.  It  is  apparent,  that  increasing  the  size  of  the 
centerline  (or  HI)  tends  to  make  the  azimuthal  mode  eigenvalues  become 
degenerate. 
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Figure  30  Resonacor  with  a hole  in  both  mirrors.  The  geometric 
output  is  through  the  annulus  of  area  it(h£  - h£)  and 
through  the  annulus  of  area  tt(A]^  - A2^) 


OB 


o 


I M-2  NO*  I 
Epi  L*3  NHI-I 

sol 


'V.  00  «|'.00  8.00  8.00  10.00 

FRF5NFI  NUMBER 


_ I 


M.00  5.00  8.00  10.00 

FRE5NFL  NUMBER 


Figure  31d  l ■ 3 


130 


NUM8CR 


FfleSWfL  MUM8FR 


Figure  32a  Eigenvalues  of  Che  resonator  of  Figure  30  for 
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In  an  attempt  to  increase  the  azimuthal  mode  discrimination  for 
KH1  - 2 and  H - 2,  the  value  of  ND  was  decreased  to  1.  Geometrically, 
this  decrease  in  D should  Increase  the  rate  at  which  the  fields  are 
coupled  over  the  outer  edge  of  mirror  2.  Since  the  higher  order 
azimuthal  mode  fields  tend  to  peak  further  out  radially  then  the  lower 
order  azimuthal  modes,  the  output  coupling  for  the  higher  order  modes 
should  Increase  more  rapidly  than  for  the  2*0  mode.  The  calculated 
eigenvalue  magnitudes  in  Figure  33  for  this  case  are  smaller  than 
those  in  Figure  32,  indicating  there  is  indeed  more  output  coupling. 
However,  there  is  no  significant  Increase  in  the  azimuthal  mode  dlserimr- 
i nation. 

The  eigenvalues  were  also  calculated  for  larger  ratios  of  HI 
and  D to  Al.  For  these  cases  the  plots  for  the  four  lowest  order 
azimuthal  modes  were  very  similar  with  the  greatest  differences  being 
between  the  2-0  and  2-3  plots.  Since  the  2-1  and  2-2  plots  lie  between 
the  2-0  and  2-3  plots,  they  are  not  shown. 

Figure  34  shows  the  eigenvalues  for  HI  - Al/2  and  D - Al/2  at  a 
magnification  M -3.  Since  D - HI  this  resonator  is  like  that  of  a 
conventional  unstable  resonator  except  the  centerline  II.  j been  expanded 
to  a radius  of  HI. 

Figure  35  shows  the  eigenvalues  for  HI  - AIM,  D-  Al/2,  and  M - 2. 
Figure  36  shows  the  eigenvalues  for  the  same  conditions  except  D - 3A1/4 
which  increases  the  area  of  the  smaller  annulus  relative  to  the  outer 
annulus.  There  is  very  little  difference  in  the  eigenvalues  for  the 
two  cases. 
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Figure  36a  Eigenvalues  of  the  resonator  of  Figure  30  for  HI  - Al/4, 
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Figure  37  shows  the  eigenvalues  for  HI  ■ Al/2,  D « 3/4A1  and 
M ■ 2.  Figure  38  Is  for  the  sane  conditions  except  D “ Al.  For 
this  last  case  where  D ■ Al  the  resonator  Is  like  that  of  resonator 
1 of  Figure  1 except  the  centerline  Is  of  radius  Al/2. 

Coopering  all  the  figures  In  this  section  shows  Increasing  the 
hole  sizes  In  the  mirrors  decreases  the  eslmuthal  mode  discrimination, 
but  Increases  the  radial  mode  discrimination.  This  suggests  an 
interesting  relationship  between  the  radial  and  azimuthal  modes  and 
also  suggests  the  following.  Since  cutting  a section  out  of  the 
resonator  along  a line  where  the  radial  coordinate  Is  a constant  Increases 
the  radial  node  discrimination,  then  cutting  sections  out  along  lines 
where  azimuthal  coordinate  Is  a constant  should  lead  to  an  Increase  In 
azimuthal  mode  discrimination.  A study  of  "wedge-coupled"  resonators 
should  be  made  to  determine  behavior  of  their  elgenmodes. 
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Figure  37a  Eigenvalues  of  the  resonator  of  Figure  30  for  HI  - Al/2, 
D - 3A1/4,  M - 2,  i - 0 


147 


SECTION  VI 


TASK  5 RESULTS 

The  purpose  of  Task  5 was  to  Investigate  the  possibility  of 
using  conformal  mapping  techniques  to  transform  circular  geometries 
to  rectangular  geometries  where  the  resonator  eigenmode  calculations 
could  be  expedited  using  the  fast  Fourier  transform.  The  conformal 
transformation  [5] 

W - exp  (Z) 

where  W - u + iv  - r exp  (j$) 
and  Z ■ x + ly 

maps,  rings  In  the  complex  W plane  into  rectangles  in  the  Z plane  as 
shown  In  Figure  39.  It  was  hoped  that  the  transform  could  be  used  to 
link  the  eigenmodes  of  strip  resonator  calculations  to  the  elgenmodes 
of  circular  resonators.  A problem  arose  which  aborted  the  exercise. 

The  transform  dictates 

r - exp(x) 

In  cylindrical  coordinates  the  field  on  mirror  2 at  r^  for  example 
depends  In  a complex  manner  on  the  field  at  each  source  point  r^.  In 
Cartesian  coordinates  the  field  at  depends  on  the  field  at  each  source 
point  x^  but  through  a Kernel  which  Is  a function  of  the  difference  in 
the  coordinates  and  x^.  Therefore  a strip  calculation  can  be  made  for 
a fixed  width  strip  resonator  anywhere  in  the  Z plane.  For  example, 
if  the  strip  mirrors  are  translated  a distance  b in  the  x direction 
relative  to  the  origin,  then  the  circular  mirrors  are  scaled  by  a factor 
exp(b)  in  the  W plane.  The  eigenmode  losses  for  the  resonator  whose 
mirrors  are  scaled  by  a factor  exp (b)  In  the  W plane  are  not  the  same  as 
the  losses  for  the  unsealed  calculation,  assuming  the  distance  between 
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the  mirrors  is  not  changed.  A possible  solution  would  be  to  scale  the 
distance  between  the  mirrors  and  the  wavelength  by  the  same  factor 
exp(b)  which  would  make  the  Fresnel  number  of  the  circular  resonator 
invarient  to  transverse  shifts  of  the  strip  resonator  as  should  be  the 
case.  The  point  has  not  yet  been  resolved. 


SECTION  VII 


SUMMARY,  CONCLUSION,  AND  FUTURE  WORK 

A computer  program  was  developed  during  Task.  1 for  calculating 
the  radial  eigenmodes  of  circular  cylindrical  resonators  where  the 
azimuthal  variation  Is  of  the  form  exp  (+  jlty)  where  l Is  any  integer. 

The  program  was  used  to  analyze  several  confocal  hole-coupled  resonator 
configurations  to  determine  whether  or  not  the  £“0  mode  vdtild  be  the 
dominant  or  lowest  loss  mode  in  the  unload  resonator.  No  purely  hole- 
coupled  resonator  was  found  in  which  the  £.*»0  mode  was  clearly  dominant 
except  for  the  uninteresting  limit  of  small  Fresnel  numbers.  One 
resonator  (Figure  lc)  was  found  in  which  all  azimuthal  modes  had 
approximately  equal  losses  (i.e. , the  modes  were  degenerate).  Loaded 
cavity  analysis  should  be  performed  using  a reasonable  gain  saturation 
profile  to  determine  the  "saturated"  mode  profile.  If  these  results 
are  promising,  an  experimental  study  of  the  resonator  should  be  under- 
taken. 

Also,  a computer  program  has  been  developed  which  demonstrates  the 
feasibility  of  a technique  for  calculating  mirror  figures  to  achieve 
specified  radial  eigenmode  profiles.  The  most  serious  limitation  is 
that  the  program  considers  only  the  radial  mode  profile,  where  in  most 
cases  of  Interest,  it  is  the  azimuthal  phase  profile  which  is  of 
paramount  interest.  The  program  should  be  generalized  to  allow  arbitrary 
radial  and  azimuthal  field  variation  so  that  the  radial  and  azimuthal 
mirror  figure  functions  can  be  calculated  to  achieve  hole-coupled  output 
with  minimal  azimuthal  field  phase  variation. 

Analysis  of  resonators  with  holes  in  both  mirrors  has  shown  that, 
when  both  the  holes  are  much  larger  than  a Frensel  number,  there  is 
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virtually  no  difference  In  the  losses  for  the  lowest  order  azimuthal 
eigenmodes,  but  that  the  radial  mode  discrimination  Is  excellent. 

Since  Increased  radial  mode  discrimination  is  achieved  by  cutting 
out  sections  of  the  resonator  along  line  of  constant  radial  coordinate 
(holes) , It  is  likely  that  increased  azimuthal  mode  discrimination 
can  be  achieved  by  cutting  out  sections  of  the  resonator  along  lines 
of  constant  azimuthal  coordinate  (wedges) . Therefore  a study  should 
be  made  of  the  eigenmodes  of  "wedge-coupled"  resonators  with  higher 
order  phase  figures. 
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APPENDIX  A 

CALCULATION  OF  THE  RESONATOR  EIGENMODES 

A computer  program  has  been  developed  to  calculate  the  bare  cavity 
scalar  radial  elgenmodes  for  exp  (j!$)  azimuthal  variation  In  circular 
cylindrical  resonators.  The  resonator  mirrors  can  have  arbitrary  radial 
curvature  within  the  Fresnel  approximation  and  arbitrary  reflectivity. 

The  elgenmodes  are  computed  numerically  via  the  Prony  method  as  described 
previously  by  Slegman  and  Miller  [1].  Their  method  Is  modified  to 
allow  for  arbitrary  mirror  curvatures,  and  the  procedure  Is  Iterated 
to  assure  convergence  and  accuracy  of  the  computed  elgenmodes  and  eigen- 
values. The  following  three  subsections  extracted  from  the  user's 
manual  for  the  program  discuss  the  relevant  mathematical  equations  and 
some  of  the  programming  details. 

1.  HAVE  PROPAGATION 

The  basic  resonator  Integral  equations  for  the  radial  field  functions 


- the  subscripts  1 and  2 refer  to  mirrors  1 and  2,  respectively 

R^  - the  phasor  radial  field  on  mirror  1,  where  the  azimuthal 

field  dependence  Is  assumed  to  vary  as  exp(j 
r - radial  coordinate 

k - wavenumber  “ 2ir/X 

X - wavelength 

B - mirror  separation 

- Bessel  function  of  order  f 

a - outer  mirror  radius 

g - resonator  psrameter  ■ 1 — B/R 

R - radius  of  curvature  • 

- single  pass  eigenvalue 
The  complete  field  Is  written  as 

- RnJlCr)  exp(j£$)  (A4) 

where  n and  1 are  the  radial  and  azimuthal  mode  Indices,  respectively. 

The  fields  (elgenmodes)  are  orthogonal  on  the  mirror  surfaces  and  can 
be  normalized  such  that 


2ir  a 

/ / E_.  (r)EnJCrir  drd<f> 


« s.0 

mn  ki 


CA5) 


For  computational  convenience  the  Integral  equations  will  be  written  In 
a different  form.  Let  us  define  a constant 


C - jl+1  ± 

J B 

a spherical  phase  factor. 


(A6) 


(A7) 


1S8 


■ 


and  a "reflection  coefficient" 


>i(rt>  - •*» 


The  letter  la  Interpreted  aa  a reflection  coefficient  Cor  function) 

2 

since  ri  /2R^  is  the  approximate  distance  the  mirror  surface  is  dis- 

2 

placed  from  planer  and  thus  kr^  /R^  Is  the  total  phase  shift  the  wave 
experiences  at  r^  upon  reflection.  For  higher  order  resonators  the  re- 
flection coefficient  can  be  generalised  to 


PjCrj)  - exp  (j 2kdtCr1)) 


(arbitrary  magnitude  variations  could  also  be  Included)  where  d^(r^)  is 
the  distance  the  mirror  Is  displaced  from  planar  (the  mirror  figure 
function) . In  terms  of  the  new  variables , the  kernel  can  be  rewritten  as 


rl»r2^  " B ) /Pj (r  j) p2 (r2) tj^ (r^ fc2 rirj 


(A10) 


Then  with  the  additional  definition 

(1)  *»1><ri> 

(r.)  - * - - — 

/Pl<rl)  /c1(r1) 


(All) 


and  the  grouping  T^Cr^)  - Pi(r1)t1(r1)t1  the  resonator  Integral  equa- 

tions may  be  written  as 


'»a\aV  ■ l2  h(r$ 

?H2)^  - ? h&r1) 


(A12) 


(A13) 


or  correspondingly  in  matrix  form  aa 


Y1S1  “ J T2S2 


(A14) 


Y2S2  - J T,S 


1“1 


CA15) 


where  l is  understood.  From  (All)  it  is  seen  that  sj^CXj)  i*  *i«ld 
on  a spherical  surface  of  radius  B Just  incident  mirror  i because  mul- 
tiplying By 

0^)  <*“> 


- exp 


transforms  the  field  to  a plane  in  front  of  the  mirror  and  multiplying 
by  transforms  the  field  to  the  mirror  surface.  Therefore,  the 

product  transforms  the  field  from  the  spherical  surface, 

to  the  mirror  surface,  and  back  to  the  spherical  surface  but  propagating 
in  the  opposite  direction. 

The  radial  fields  on  the  spherical  surfaces,  can  be  normalized  such 

that 


/ .<«  (r)  s Cr)  T(r)  dr  - (A17) 

where  the  Integral  is  over  the  mirror  surface. 

Finally  the  composite  matrix  equation  for,  say  the  fields  on  mirror 
1,  may  be  written 

yS1  - J T2  J Tj.  Sx  (A18) 

where  the  round  trip  eigenvalue  y * Y2  Y2*  Yor  8Ba11-  number  of  points 
on  both  mirrors  it  is  useful  to  form  the  round  trip  operator  matrix 
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M - JTj  JT^.  However,  as  the  number  of  field  grid  points  Increases,  two 
computational  savings  can  be  realized  by  not  computing  M.  First,  It  takes 
one  fourth  the  computer  storage  to  store  the  the  real  symmetric  matrix  of 
Bessel  functions  J as  opposed  to  the  complex  matrix  M.  Second,  the  nuaber 
of  multiplications  required  to  fora  M can  exceed  the  number  of  multiplica- 
tions required  to  find  the  dominant  eigenmodes  by  the  Iterated  Prony 
method  which  will  be  described  In  the  next  section. 

2.  ITERATED  FRONT  METHOD 

The  basic  Prony  calculation  used  Is  exactly  aa  described  by  Slegman 
and  Miller  [1]  except  the  procedure  Is  repeated  until  convergence  Is 
achieved  and  the  nuaber  of  eigenvalues  Is  changed  as  required  to  facil- 
itate convergence.  The  following  discussion  briefly  describes  the  com- 
putational procedure. 

Formally,  we  wish  to  find  the  first  n dominant  eigenvectors  and 
eigenvalues  1 ■ 1 to  n of  the  matrix  equation 

MU  - Uy  (A19) 

where  M Is  a complex  symmetric  matrix.  In  resonator  terminolgy  the 
eigenvectors  are  often  referred  to  as  "eigenmodes,"  "modes,”  or  "fields." 
An  Initial  vector,  or  field  distribution,  Vq  Is  constructed  and  propagated 
n times  through  the  resonator  using  the  wave  propagation  procedure  de- 
scribed In  the  previous  section.  The  higher  order  modes,  smaller  eigen- 
values, are  attenuated  more  rapidly  than  the  lower  order  modes,  which 
have  larger  eigenvalues.  The  Initial  field  and  sequence  of  propagated 
fields  may  be  written  as: 
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Vo  - Vl  + V2  + k3°3  + 


V1  " " Wl  + V2°2  + V3U3  + 


Vn  " "Vl  * VX  + VS°2  + *3^3  + 


(A20) 


where  initially  the  k^'s  and  y^'s  are  unknown.  Two  matrices  of  dot 
products  are  generated 


(A21) 


(A22) 


where  r » n - 1 and  the  dot  products  are 


- Z Cv.).  Cv.) 


I k ^j'k 


(A23) 


The  F matrix  Is  Inverted  and  the  eigenvalues  are  computed  as  the  roots 


of  the  polynomial 


Y°  + QrYn_1  + • • - QjY  + Q0  - 0 


(A24) 
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Using  the  k^s,  y^s,  end  V^s,  i - 1 to  n,  estimates  U^s  of  the 
U^s  are  computed  from  CA20) . Since  there  may  have  been  components  of 
higher  order  modes  than  n In  the  original  vector  Vq,  we  do  not  at  this 
point  know  how  accurate  our  computations  of  the  elgenmodes  are.  There- 
fore a new  Initial  vector  Is  formed,  using  the  just  computed  eigenmode 
estimates  as 


V - U.  + U,  + . . . U 
o l L n 


(A30) 


where  the  are  normalized  Such  that 


l ■ 1 

k 


CA31) 


The  procedure  Is  repeated  until  the  calculated  expansion  coefficients 
k's  are  unity  and  the  eigenvalues  do  not  change  from  one  calculation  to 
the  next.  Ideally  if 


frn+lJ  " frj 


(A3  2} 


the  procedure  will  always  converge,  hut  when 


fo+ll  * Ki I 


(A33) 


convergence  can  take  a large  number  of  Iterations.  When  this  situation 
occurs,  it  Is  expedient  to  change  n.  If  the  eigenvalues  tend  to  separate 


as  n Increases,  as  Is  the  case  in  resonator  calculations,  n should  be  in- 
creased. However  n cannot  always  be  increased  since  the  finite  word 
length  of  the  computer  will  eventually  truncate  components  of  the  higher 
order  mode  in  Vq  in  equation  (A20) , i.e., 

Gr) 

will  become  smaller  than  the  resolution  limit  of  the  computer. 

3.  FAR-FIELD  CALCULATIONS 

By  convention  in  the  computer  program,  the  output  end  of  the  res- 
onator is  at  mirror  2.  For  display  purposes,  the  field  incident  on 
mirror  2 is  transformed  to  a plane  via 

u(*2)  * /t(r2)s(r2)  (A35) 

The  far-fleld  pattern,  f(r).  Is  computed  by  taking  that  part  of  u(r2> 
which  exits  the  resonator,  ut^r2),  and  propagating  it  to  the  focus  of 
a lease.  For  convenience  the  focus  of  the  lense  is  B,  the  cavity  length, 
so  the  previously  computed  Bessel  function  can  be  used.  The  field  at 
the  focus  of  the  lens  is  computed  from 

fCr)  - C / ut^r2^r2dr2  (A36) 

where  the  integration  is  over  the  output  plane.  The  integrated  far- 
field  Intensity,  p(r),  can  then  be  computed  as 

p(r)  - £r  IfCr'JlVdr*  (A37) 
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APPENDIX  B 


LEASTS  SQUARES  METHOD  FOR  CALCULATING  RESONATOR 
PARAMETERS  TO  ACHIEVE  A DESIRED  EIGENMODE  PROFILE 

A non-linear  regression  algorithm  has  been  utilized  to  calculate 
radial  mirror  figures  of  open  resonators  such,  that  the  radial  profile 
of  the  lowest  loss  eigenmode  Cfor  a particular  azimuthal  symmetry)  1s 
the  best  fit  in  a least  squares  sense  to  a user  specified,  "desired", 
eigenmode  profile.  The  method  used  is  that  of  Marquardt  [7]  which  is 
a maximum  neighborhood  method  for  least  squares  estimation  of  nonlinear 
parameters.  It  performs  an  interpolation  between  the  Taylor  series 
method  and  the  gradient  method  in  an  attempt  to  combine  the  best  features 
of  the  two  methods  while  at  the  same  time  avoiding  their  most  serious 
limitations.  The  algorithm  described  here  for  ease  of  reference  is 
identical  to  that  of  reference  [8]  except  for  changes  in  nomenclature 
and  the  inclusion  of  point  error  weighting. 

The  model  to  be  adjusted  to  the  data  can  be  written  in  the  form 

U - fO:;81,82,  . . . Bk)  (Bl) 

where  r is  the  independent  variable(s)  and  8^,  82,  * * • are  the 

/v 

population  values  of  the  k parameters  to  be  estimated  and  U is  the 

A 

dependent  variable.  In  the  context  of  the  resonator  problem  U is  the 
lowest  loss  radial  eigenmode  profile  (or  intensity)  for  a particular 
axlmuthal  symmetry,  r is  the  radial  field  coordinate,  and  the  6's  are 
the  resonator  parameters  which  one  wishes  to  vary.  The  observed  points 
will  be  designated  as 

UC)^)  1 - 1,  2,  ...  N (B2) 


166 


r 


which  for  the  cases  investigated  to  date  are  the  desired  radial  field 
points  in  a plane  adjacent  to  one  of  the  resonator  mirrors  or  desired 
radial  points  in  the  Integrated  far-field  intensity  profile.  Standard 

A 

techniques  can  be  used  to  calculate  the  radial  mode  profile  U at  the 
specified  coordinates  for  a given  resonator  geometry  and  radial  mirror 
figures  which  are  dependent  on  the  g's.  The  problem  could  be  structured 
even  more  generally  so  as  to  optimize  azimuthal  variation  In  the  mirror 
figures  provided  the  total  field  was  calculated  without  restriction  on 
the  type  of  azimuthal  symmetry. 

The  problem  is  to  compute  those  estimates  of  the  parameters  which 
will  minimize 


$ - E W.  [u.  - u.]‘ 
i-1  1 1 


(B3) 


i -3 

. 

1 


where  u^  is  the  value  of  uCr^)  predicted  at  the  i-th  data  point  and 
is  the  associated  weight  for  that  point  error,  (u-magnitude  and/or  phase  of  U) 
The  Taylor  series  aspect  of  the  combined  Marquardt  Method  will  now 
be  described.  Writing  the  perturbed  model  in  a Taylor  series  through 
the  linear  terms  gives 

k 


< utoj,  f (rt,  y + £ ||l  (3t)J 


(B4) 


or 


< u > - f + P6 
-o  — t 


(B5) 


In  (B4) , £ is  replaced  by  b^  the  converged  value  of  b^  being  the  least 
squares  estimate  of  8 and  f±  is  the  value  of  the  function  evaluated  at 
the  i-th.  data  point.  The  vector  5^  is  a small  correction  to  b_,  where 
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th«  subacripC  C is  used  Co  designate  Che  £ was  calculated  using  Che 
Taylor  series  method,  and  P is  an  [n  x It]  matrix  whose  elements  are 
3ft/3bj.  The  brackets  <>  are  used  to  distinguish  predictions  based  on 
the  linearized  model  from  those  based  on  the  actual  non— linear  model* 
Therefore,  the  estimated  sum  of  squares  error  CSSE)  is  designated  as 


<♦>  - E W. [u.  - <u  >]  ( 

i-1  1 1 

In  order  to  minimize  «|>>  the  critical  point Cs)  must  be  found,  l.e.. 


“ 0,  j - 1,  2,  ....  k 


which  requires  from  (B4)  and  (B6)  that 

<*>  - 2 Vui  - *i  "2  f£*  C«t)J2 

i-1  1 1 1 j-1  3bj  t 3 


Z W. [u  - f.  - Z - 0 

i-1  1 L 1 j-1  3bj  C 3 3bi 


n k 


ifl  ” fi*  3^ 


E Z W 5ft  Cfi  ) 
i-1  j-1  1 9bj  36i  C 3 


or  in  matrix  notation 


A *1-M. 


(Bll) 


where 


AIt  x k]  . pTr 


(B12) 


[n  x k]  _ 


||i).  1-1.2. 


.,  n;  j - 1,  2,  . . .,  k (B13) 


l‘k  1 11  ■ j,  V“l  - 'J  If}  • J - 1.  2.  • 


. .,  k (B14) 


- PaD(u  - f) 


(B15) 


and  t Is  Che  new  perturbation  as  Influenced  by  the  weighting  and 
D[n  x ®1  ^ a diagonal  matrix  of  the  weights.  In  practice  is  only 
corrected  by  a factor  of  6^  so  that  extrapolation  does  not  go  beyond 
the  region  where  f Is  adequately  described  by  CB4) . Failure  to  do  so 
would  cause  divergence  of  the  iterates. 

The  steepest  descent  method  takes  a step  In  the  direction  of  the 
negative  gradlaent  of  4 from  the  current  trail  value.  Thus, 


**  1b2 3b  J 


CB16) 


Here  again  the  step  size  requires  careful  control  to  insure  convergence. 

Having  noted  the  Inadequacies  of  both  the  Taylor  and  the  steepest 
descent  methods,  note  also  that  any  proper  method  must  result  in  a cor- 
rection vector  that  is  within  90°  of  the  negative  gradient  or  else  $ 
would  Increase  with  the  iterate.  Marquardt’s  Method  Is  an  interpolation 
between  5^  and  <5^  for  monotonic  convergence  of  the  procedure.  This 
method  has  its  foundation  in  three  theorems  which  will  be  stated  without 
proof.  Proof  of  these  theorems  may  be  found  in  Marquardt  [7]  or  Meeter 
[8].  Meeter  states  and  proves  a somewhat  stronger  version  of  Theorem  1. 

Theorem  1:  Let  X>0  be  arbitrary  and  let  6 satisfy  the 

equation 


(A  + XI)6  - i 


(B17) 
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then  5 minimises  <♦>  on  the  sphere  whose 
radius  1 1 5 1 1 satisfies 

inn2-  mu2 . (Bi8> 

Theorem  2:  Let  5(X)  be  the  solution  of  (B17)  for  a given 
value  of  X.  Then  ||6,(X)||*  is  a continuous 
decreasing  function  of  X,  such  that  as  X-*», 
l|5(X)||a-0. 

Theorem  3:  Let  y he  the  angle  between  6^  and  5 . Then  y 

is  a continuous  monotone  decreasing® function  of 
X such  that  as  X-***,  y-K) . Since  5 is  Indepen- 
dent of  X,  It  follows  that  5 rotates  toward 
6 as  X-*-.  ~° 

Since  the  Marquardt  Method  combines  both  the  Taylor  series  method 
and  the  steepest  descent  method.  It  Is  necessary  to  note  one  relevant 
property  of  these  methods . The  solution  of  (B11J  is  invariant  under 
linear  transformation  of  the  h-epace.  However,  the  steepest  descent 
solution  is  not  scale  invariant.  It  becomes  necessary  to  scale  the  b— 
space  in  some  convenient  manner.  Therefore,  the  h-space  is  scaled  in  units 
of  the  standard  deviations  of  the  derivatives  Sf^Sb^,  taken  over  the 
sample  points  1 “■  1,  2,  . . .,  n.  This  choice  of  scale  causes  the  A 
Batrlx  to  be  transformed  into  the  matrix  of  simple  correlation  coefficients 
among  the  ^fj/Sb^. 

Thus,  the  scaled  matrix  A*  and  the  scaled  vector  j*  become 
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and  the  Taylor  aarlea  solution  leads  to 


A fi  m a* 
— t 


(B21) 


Then 


si  - «!'  ^ 


(B22) 


The  algorithm  la  now  clear.  Specifically,  at  the  r-th  Iter- 
ation the  equation 

(A*(t>  + X<r)  I)«*(r)  - * (r> 


(B23) 


la  cona  true  ted.  Thla  equation  la  aolved  for  and  (B22)  la  uaed 


t(r) 


to  obtain  & ' 1 . The  new  trial  paraaieter  vector 


b(r+!)  . b(r)  + fi(r) 


(B24) 

fr+1) 

will  lead  to  a new  sum  of  aquares  error  4V  . It  la  eaaentlal  to 
(r) 


select  \ such  that 

< *<r>  (g25) 

It  la  clear  from  the  foregoing  theory  that  a sufficiently  large 

always  can  be  found  such  that  (B25)  la  satisfied,  except  when 


.(r) 


b'  is  already  a minimum  of  ♦.  Therefore,  a trial  and  error  procedure 


(r) 


la  required  to  find  a value  of  X'  such  that  (B25)  la  satisfied  and 
rapid  convergence  Insured,  l.a.  minimize  * in  the  (approximately) 
maximum  neighborhood  over  which  the  linearised  function  will  give 
adequate  representation  of  the  non-linear  function. 

The  strategy  employed  is  as  follows: 

Let  v > 1. 


Let  X^  denote  the  value  of  X from  the  previous  itsratlon 


Initially,  let  X ^ - 10~2,  for  example. 

Compute  *(X^r“^)  and  *(X^r_1Vv)  (see  note  below). 


m 
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i.  If  *(X<t“1)/v)  < *<r),  Ut  X(r)  • X^Vv 
li.  If  *(X(r“1)/v)  > *(r),  and  *(X(r"1))  <*(r),  lat 
X<r>  - X^15. 

111.  If  ♦(X^r”1fv)  > and  • (X^r*”^)  > Increase  X 


by  successive  multiplication  by  v until  for 
«.  ♦a<t-l)v“)  < «(r).  ut  xM  - 
NOTE:  If  X^-1)  1.  already  negligible  by  coaparlaon  vlch  1.0  to 


smallest 


the  number  of  significant  figures  carried,  then  go  to  test 
11  or  111  Immediately  without  computing  ♦(A^r~^/v) , and 
Ignore  comparisons  Involving  *(X^r”^/v). 


On  occasion  In  problems  where  the  correlations  among  the  parameter 


estimates  are  extremely  high  ( > 0.99),  It  can  happen  that  X will  In- 


crease to  unreasonably  large  values.  It  has  been  found  helpful  for 


these  Instances  to  alter  test  111.  The  revised  test  Is: 


J,<r41>  . b(r>  + K(r)£(r)>  K<r)  < , 


Noting  that  the  angle  la  a decreasing  function  of  X^r\  select  a 


criterion  agnle  YQ  < j and  take 


:(r)  - 1 if  y(r)  > y 

— ' C 


However,  If  test  111  Is  not  psssed  even  though  ' has  been  Increased 


until  Y^  < YQ»  then  do  not  Increase  X^  further,  but  take 


sufficiently  small  so  that  < 4^,  This  can  always  be  done 

since  y(r)  < YQ  < f [71. 

The  procedure  is  considered  converged  when( | 6 ^ ^ | /(t+|bj ^ | ) < e, 
for  all  j ■ 1,  2,  ...,k  and  suitable  choices  for  e and  t,  for  example 


•5  _3 

10  and  10  , respectively. 
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